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Abstract 

The interacting boson model is extended to describe the spectroscopy of 
superdeformed bands. Microscopic structure of the model in the second mini- 
mum is discussed and superdeformed bosons are introduced as the new build- 
ing blocks. Solutions of a quadrupole Hamiltonian are implemented through 
the 1/N expansion method. Effects of the quadrupole parameters on dynamic 
moment of inertia and electric quadrupole transition rates are discussed and 
the results are used in a description of superdeformed bands in the Hg-Pb 

and Gd-Dy regions. 
21.60Fw 
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I. INTRODUCTION 



Spectroscopy in the super deformed minimum has reached a certain level of maturity 
to justify a phenomenological analysis of the available data (see for recent reviews). 
Such an approach would be useful in systematizing the data and would also provide a 
complimentary perspective to the more microscopic theories. For this purpose, we use the 
interacting boson model (IBM) || which has been established as one of the simplest and 
most versatile collective models. It has been especially successful in correlating spectroscopic 
data in deformed nuclei in terms of a few parameters of a quadrupole Hamiltonian [||. 

Microscopic study of the nucleon pair structure in the superdeformed well || indicates 
that, compared to the deformed nuclei, they have about three times more active pairs 
of nucleons (bosons), and the L = 4 pair (g boson) plays a much more significant role. 
As numerical diagonalization of an sdg-lBM Hamiltonian is not possible for more than 10 
bosons, one needs alternative methods of solution to apply the IBM to superdeformed nuclei. 
Here, we use the angular momentum projected mean field theory which leads to a 1/N 
expansion for all matrix elements of interest |J. Accurate representation of high-spin states 
in the 1/N expansion formalism requires terms to order (L/N) 6 which have been obtained 
recently using computer algebra [|7j . The extended formalism provides an analytical method 
for analysis of superdeformed states which is both accurate and efficient. 

The plan of the paper is as follows. After reviewing the microscopic basis of the IBM 
for superdeformed nuclei, we introduce the 1/N expansion formalism and discuss its recent 
extension to higher orders. We then use the 1/N expansion formulas for a quadrupole 
Hamiltonian to study systematic features of dynamic moment of inertia and B(E2) values. 
The results are used in a description the superdeformed bands in the Hg-Pb and Gd-Dy 
regions. 
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II. MICROSCOPIC BASIS 



In this section, we study the typical structure of strongly deformed states and investigate 
the relation between the superdeformation and the IBM based on ||. For this purpose, we 
use the Nilsson + BCS model with particle number projection. Superdeformed states can be 
characterized as ground states in a superdeformed potential well which is separated from the 
normal one by a potential barrier. For such ground-like states which show strongly collective 
nature, this model seems to work well. Using the experimental deformation parameters 
and electric transition probabilities (or moments) as input, one can obtain reasonable wave 
functions. These wave functions are analyzed from the viewpoint of collective nucleon pairs, 
which leads to a natural extension of the usual IBM. 

We briefly summarize the formulation of the Nilsson + particle-number-conserving BCS 
model. The single particle orbits in a deformed potential are described well by the Nilsson 
Hamiltonian || 

Tl^ 771 4 

^Nilsson = -T^V 2 + —LUy[l ~ -5P 2 (C0S6)} 

— 2KuoqkI ■ s — hujQKfi(l 2 — (1 2 )n)i (1) 

where 5 is the deformation parameter and P2(cos8) denotes the Legendre polynomial. The 
term (I 2 ) jy = ^N(N + 3) is the expectation value of I 2 averaged over one major shell with the 
principal quantum number N = 2n + I. The value of the oscillator frequency for a mass-A 
nucleus is determined from Tiujq = 41A~^, which is about 7.1 MeV for the superdeformed 
nuclei in the Hg-Pb region (A ~ 190) and 7.7MeV for those in the Dy-Gd region (A ~150). 
We use the usual values for the parameters k and fi which are 0.0637, 0.60 for proton 
orbits and 0.0637, 0.42 for neutron orbits, respectively pj. In order to include short range 



correlations, the monopole pairing interaction is added to the Nilsson hamiltonian |L0 

H = ^Nilsson + GPtP, (2) 

where G denotes the pairing strength parameter, and 
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^ t = E 44 . (3) 

fe>0 

is a pair creation operator. Here c\ stands for the creation operator of a nucleon in the 
spherical single particle orbit k, and k denotes the time reversed state of k. This Hamiltonian 
is solved by the variation using a BCS wave function 

I *) = n( u « + u «°a4) I °>> ( 4 ) 

a>0 

where denotes the creation operator for a nucleon in the deformed canonical (Nilsson) 
orbit labeled by a. The particle number conservation has been found to be important in 
the case of weak pairing correlations and also for moments of inertia of high spin states in 



the cranking calculation of the superdeformed states [|TT|J12[] . Thus we carry out the particle 



number projection before variation according to the method given in Ref. [ 13| . The solution 
corresponds to the minimum of the number projected energy 
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E F m = h ', ' / , (5) 



where P N denotes the particle number projection operator. 

The deformation parameters of the superdeformed states in the Hg-Pb (Dy-Gd) region 
are given by 5 ~ 0.40 (0.50), which is equivalent to the axis ratio of 5:3 (2:1). Because of 
this strong deformation, it is insufficient to take only one active major shell and take into 
account the corrections due to the core-polarization effect through renormalization in one 
major shell. Thus we first seek a suitable model space for description of superdeformed 
states. For simplicity, we turn off the pairing force which is not important for this purpose. 
We take 194 Hg (JV=114, Z=80) and 152 Dy (iV=86, Z=66) as examples of the Hg-Pb and 
the Dy-Gd regions, respectively. 

In order to define the model space necessary for description of superdeformed states, we 
utilize the intrinsic quadrupole moment Qq. The intrinsic quadrupole moment is calculated 
in the space of all spherical orbits up to the principal quantum number N = iV max . Then 
-^Vmax is increased until the value of Q is saturated to a good extent. From this procedure 
we obtain iV max = 12. The corresponding values of Qo for proton and neutron orbits are 



196 and 296 for 194 Hg, and 196 and 276 for 152 Dy, respectively. Note that the experimental 
values of Qq are 18 ± 3 eb for 152 Dy |14j] and 17 ± 2 eb for 194 Hg [[OJ, which are consistent 



with the present results if we take the bare charges, e p = 1 and e n = 0. 

We now consider the inert core of super deformed states. The Nilsson wave function is 
obtained by putting all nucleons in the Nilsson orbits from the bottom. One Nilsson orbit 
can be expanded as a linear combination of many spherical harmonic oscillator orbits, and 
the square of expansion coefficients gives the occupation probability of each spherical orbit. 
We expand all the occupied Nilsson orbits and sum up all the occupation probabilities which 
belong to the same spherical orbits, to obtain the total occupation probability for a given 
spherical harmonic oscillator basis. Due to the strong quadrupole field, one Nilsson orbit 
spreads over many spherical orbits. Thus the orbits with very high single particle energy 
can gain some finite occupation probabilities, while the occupation of the orbits with small 
single particle energy may become incomplete. Nevertheless several lower spherical orbits 
are occupied almost completely and can be considered as a new inert spherical core for the 
superdeformed states. Note that we do not take the usual "hole" picture as it is meaningful 
only for states whose configuration are well described within one major shell. 

First consider the case of 194 Hg. In Fig. |I] the occupation probability of each spherical 
harmonic oscillator orbit is shown for neutrons (a) and protons (b). The orbits are ordered 
according to their single particle energy at 5 = as lsi/ 2 , IP3/2, • • •• The case of 5 = —0.13 
which simulates the deformation of normal oblate states is also shown for comparison. For 
normal deformation, it is seen from Fig. p]-a that the occupation of the proton orbits is 
almost complete at 2d 5 / 2 (Z = 64), while the occupation probability is almost vanished for 
orbits above Z = 82. These results suggest that we can consider the Z = 64 subshell as an 
inert core and three valence orbits (2d 3 / 2 , 3si/2 and l/in/2) as active. This gives the valence 
proton number as Z v = 16. In the case of 5 = 0.40, the proton orbits are almost completely 
occupied up to Z = 50. Above Z = 50, the occupation probability drops suddenly though 
it remains about 10% over many orbits. Clearly, one should incorporate the contributions 
of these high energy orbits. Thus it is reasonable to take the Z = 50 spherical inert core 



and include quite many orbits above there as active valence orbits. In this case the valence 
proton number becomes Z v = 30. 

In the same way, it can be seen form Fig. |I|-b that the spherical inert core for neutron 
orbits are N = 100 (N v = 14) and 82 (N v = 32) for 5 = -0.13 and 0.40, respectively. The 
active valence orbits for 5 = —0.13 are 2/ 5 / 2 , 3p 3 / 2 , 3pi/ 2 and U13/2, while for 5 = 0.40 it is 
still insufficient to include only two or three major shells. 

We can see a similar behaviour of occupation probabilities in the wavefunctions of 152 Dy, 
which is shown in Figs. |l|-c for proton and |l|-d for neutron orbits. In these figures two cases 
of 5 = 0.50 (super deformed state) and 0.25 (normal prolate state) are compared. It is clear 
that Z = 50 and N = 82 inert cores are good for normal states, while Z = 40 and N = 50 
cores are suitable for superdeformed states. Thus the inert core of superdeformed states 
becomes much smaller than that of normal states in both the Hg-Pb and Dy-Gd regions. 

From the viewpoint of the IBM, the number of bosons is determined by half of the 
number of valence nucleons. Because of the small inert core the number of bosons increases 
significantly for superdeformed states in comparison with that in the usual IBM. In fact in 
the case of 194 Hg, the boson number in the usual IBM is N noima i = (82 — 80)/2 + (128 — 
114) /2 = 8 by taking the usual hole picture, and iV normal = (80-64)/2 + (114- 100)/2 = 15 
with the particle picture mentioned above. On the other hand, the number of bosons for 
super deformation becomes A^ super = (80 — 50)/2 + (114 — 82)/2 = 31. Similarly, in the case of 
152 Dy, iV normal = (66-50)/2 + (86-82)/2 = 10 while iV supcr = (66-40)/2 + (86-50)/2 = 31. 
In general, iV SU p er is about three times larger than iV norma i. Note that the number of proton 
bosons and neutron bosons are close in these two cases, and this approximate equality 
seems to be a general tendency of the superdeformed states. This result can be naturally 
understood since equal numbers of valence protons and neutrons maximizes the attractive 
proton-neutron interaction. 

Next we consider the effects of pairing correlations on the structure of wave functions of 
superdeformed states. The strength parameter G of the pairing interaction should be chosen 
depending on the model space. Since the value of G for such a large space is not known 
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empirically, we first describe normal states within the extended valence space, and determine 
the value of G by requiring that the pairing gap A takes a reasonable value. For 194 Hg the 
value of G has turned out to be 0.06MeV which gives A ~ 1 MeV. Using this value, we 
investigate the effect of the pairing correlations on the structure of valence wave functions. 
For this value of G, the gap for superdeformed states becomes about A = 0.5 MeV for both 
proton and neutron orbits. In contrast to normal deformed states, which are sensitive to 
changes in values of G, the superdeformed states are almost insensitive to G values (the 
intrinsic quadrupole moment and the occupation probabilities change very little). Thus the 
following discussion about the structure of valence wave function of super deformation is 
almost independent of pairing correlations. 

We can investigate the relation between the superdeformation and the IBM by analyzing 
valence wave functions from the viewpoint of collective nucleon pairs. Since the bosons in 
the IBM are understood as images of these pairs, such an analysis is essential in establishing 
a microscopic basis for the super IBM. We consider 194 Hg as an example. 

The Nilsson + particle-number-conserving BCS wave function can be expressed as the 
condensed state of coherent Cooper-pairs in the deformed potential (TB| 



*>« (4)^(4)^10), (6) 



acting on the inert core |0). In this expression, 4 (4) denotes the creation operator of 
a Cooper-pair in proton (neutron) orbits and N n (N v ) means half of the valence proton 
(neutron) number. These yl-pairs can be decomposed into a linear combination of collective 
nucleon pairs with good angular momenta 

yi f = XoS^ + x 2 D\ + x A G\ + ■■■, (7) 

where S\ D>, G', ■ ■ ■ denote the collective nucleon pairs with spin-parity J 71 = + , 2 + , 4 + , 
• ■ • and the xj's are amplitudes. The probability of each pair in the /t-pair is given by the 
square of each amplitude, and is listed in Table | for two cases of S = —0.13 and S = 0.40. 
It is well known that in the case of normal deformation the dominant components are the 



S- and D-pairs [17 , 18| . In fact, these two components account for 100% probability in the 
case of 5 = —0.13. In the case of 5 = 0.40, the total probability of the S- and -D-pairs is 
about 80% and we can conclude that these pairs are still dominant in the /1-pair. However 
the probability of the G-pair is now sizable, and it can no longer be neglected in a detailed 
description of high-spin states. It should be noted that the ratio of the S'-pair to the other 
pairs is quite similar to that of s-boson to the other bosons in the SU(3) limit of the IBM, 
which are shown in the same table. This suggests that the SU(3) limit of the sdg-lBM could 
provide a reasonable phenomenological framework for superdeformed states. 

To summarize the microscopic results, we emphasize two important points for the de- 
scription of superdeformed bands in the IBM: One is the significant increase in the boson 
number, and the other is the importance of g-bosons. In addition, it has been found that 
the bosons for superdeformed states carry the collectivity over many major-shells and that 
the SU(3) limit is a reasonable starting point. 



III. 1/N EXPANSION FOR SUPER IBM 

A simultaneous description of the spectroscopy of normal and superdeformed states re- 
quires rather complicated wave functions, therefore we focus on the latter here and leave the 
complete picture for future work. We introduce the superbosons s, d, g as the boson images 
of the S, D, G collective nucleon pairs in the superdeformed well (bold face notation is used 
for super bosons to distinguish them from the normal ones). The quadrupole Hamiltonian 
for this system of bosons has the form 

H = -kQ ■ Q, (8) 

where the quadrupole operator is defined as 

Q = [s+d + d+sp) + g 22 [dtdp) + g 24 [dtg + gtdp> + g 44 [g t g] (2) . (9) 

Here brackets denote tensor coupling of the boson operators and b\ m = (— l) m bi_ m . The 
parameters qji in Eq. (||) determine strengths of boson interactions relative to the s — d 
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coupling. Since the SU(3) limit is used as a reference point in the rest of the paper, we 
quote the values for the quadrupole parameters in this limit; g 2 2 = —1-242, g 2 4 = 1-286, 
g44 = —1.589. As stressed in the introduction, numerical diagonalization of this Hamiltonian 
for N ~ 30 bosons is not possible even on a supercomputer. The large number of bosons are, 
however, an advantage for the analytic 1/N expansion technique which we employ here for 
solving the Hamiltonian Eq. (§). The 1/N expansion method has previously been discussed 
in detail || and the recent extensions to higher orders are given in Ref. |7fl. Therefore, 
we give only a short account of the formalism here, focusing mainly on the accuracy of 
the results for high-spin states. The starting point of the 1/N calculations is the boson 
condensate 

|iV,x) = (iV!)- 1 /2(bt)^|0), bt = x st + x 2 d\ + x 4g J, (10) 

where x\ are the mean field amplitudes to be determined by variation after projection (VAP) 
from the energy expression 

E L = (N, X \HP L \N, X )/(N,x\P L \N,x). (11) 

Here Pq denotes the projection operator. The resulting energy expression is a double 
expansion in 1/N and L = L(L + 1), and has the generic form 

— 2 £^(^)", < 12 > 

where a = J2i l x i an d the expansion coefficients e nm involve various quadratic forms of the 
mean fields X[, e.g., e 00 = (^2ji(iOK)\2Qi)qjiXjXi) 2 . The coefficients e nm have recently been 
derived up to the third order, (L/N 2 ) 3 , using computer algebra 0. 

Another observable of interest in the study of superdeformed states is the E2 transitions. 
Assuming that the quadrupole transition operator is the same as in the Hamiltonian, i.e. 
T(E2) = eQ where e is an effective boson charge, the E2 matrix elements are given by 

(L || T(E2) || L - 2) = eNL(L0 20\L - 2 0) [mi + m 2 L(L - 1) /N 2 ] (13) 
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where L = [2L + l] 1 / 2 and the coefficients m n are given in Ref. 0. The first term in Eq. 



(|13|) gives the familiar rigid-rotor result. The second term is negative and is responsible for 
the boson cutoff effect in E2 transitions. 

Before applying the 1/N expansion results, we compare them with those obtained from 
an exact diagonalization of the Hamiltonian. Diagonalization is carried out for N = 10 
which is the maximum possible boson number for this purpose. The quadrupole parameters 
?22j #24, #44 are scaled down from their SU(3) values with q = 0.7 which gives an adequate 
parametrization for the Hg-Pb region. Fig. Q-a compares exact results for the dynamic mo- 
ment of inertia J 1 ^ (circles) with three different 1/N calculations. The solid line shows the 
third order VAP results which is seen to follow the exact results very accurately. The second 
order VAP (dotted line) and the third order VBP results (dashed line) break down around 
spin L ~ 2N. Hence for description of high-spin states, the third order 1/N expansion 
formulas with VAP seem to be both necessary and sufficient. In Fig. §-b, the exact E2 
transition matrix elements (circles) are compared with those obtained from Eq. ([H|) (line). 
The agreement is again very good up to very high-spins. Note that the accuracy of the 1/N 
expansion results improves with increasing N, hence in actual applications with N ~ 30, 
one would expect an even better agreement. The test case discussed here indicates that 
the extended formalism can be applied with confidence in the spin region L = N-3N which 
covers the presently available data range for superdeformed bands. 



IV. APPLICATIONS TO SUPERDEFORMED BANDS 

In this section, we apply the 1/N expansion formulas first in a systematic study of 
dynamic moment of inertia and B(E2) values, and then to describe the experimental data 
on superdeformed bands. Since k is a scale parameter for energies, we need to study the 
effect of the three quadrupole parameters ^22, 924, 944- Fig. ^| shows the effect of variations 
in each qji on dynamic moment of inertia while the other two are held constant at the SU(3) 
values. Here q denotes the scaling parameter from the SU(3) values. Thus q = 1, corresponds 
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to the SU(3) limit which exhibits the rigid-rotor behaviour. To describe the variations in 
one needs to break the SU(3) limit. From Fig. |3] it is seen that is most sensitive 
to #24 (note the different scales in the three figures). The other (diagonal) parameters have 
smaller and opposite effect on . Since the amount of data does not justify use of too 
many parameters, we prefer to scale all three with the same parameter q. The result of this 
simultaneous scaling is shown in Fig. |]-a which is essentially the same as the one for §24 in 
Fig. H An interesting feature of these results is that the quadrupole Hamiltonian has the 
scope to describe both the increases and decreases in . For q < 1, the s — d coupling 
is relatively stronger than the d — g coupling which results in loss of monopole pairing with 
increasing spin, and hence increase in . The opposite happens for q > 1. Fig. [|-b shows 
the effect the simultaneous variations in the quadrupole parameters on the B(E2) values. 
The curving down of lines is due to boson cutoff which is most effective for smaller values 
of q. 

In the light of the above systematic studies, we have carried out fits to the available data 
on superdeformed bands in the Hg-Pb and Gd-Dy regions. The boson number is determined 
from microscopies, and k and q are fitted to the data. The parameter values are given in 
the figure captions and the data are taken from the compilation in Ref. [I5|. Figs. [5] and 



P compare the experimental dynamic moment of inertia (circles) with the calculated ones 
(lines) in Hg and Pb isotopes, respectively. In all cases exhibits a smooth increase which 
is well reproduced by the calculations. The situation in the Gd-Dy region is not as favorable 
for our simple collective model as the other region, because there are definite signs indicating 
the importance of the single particle degree of freedom. For example, in 144-146 Gd there are 
sudden jumps in which are probably due to particle alignment effects p0[ . In 148 ~ 150 Gd, 
behaves reasonably smoothly so we have attempted to describe them (see Fig. [7]). The 
average behaviour of in 148 Gd is reproduced but the model fails in the case of 150 Gd, 
underscoring the importance of single particle degree of freedom. For a better description 
of the data, one needs to incorporate particle alignment effects in the present formalism 
by including two-quasiparticle states in the basis [21]. The dynamic moment of inertia of 
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superdeformed bands in Dy isotopes exhibit an entirely different behaviour (Fig. |8]). They 
are very close to the rigid- rotor values, and hence the SU(3) limit as reflected in the values 
of q ~ 1. 

The B(E2) values provide a complimentary observable to which could be used as 
a further test of the model. In Fig. |9], the available B(E2) data in 192 ~ 194 Hg (circles) are 
compared with the calculations. A reasonable description is obtained using boson effective 
charges e = 0.12 — 0.14 eb which are typical values used in the normal IBM calculations. 
Since the B(E2) values are sensitive to the boson number (they vary as N 2 ), this provides a 
consistency check on the microscopically derived boson numbers. A further N dependence is 
provided by the boson cutoff term in Eq. (|T3|) which causes a drop in the calculated B(E2) 
values at high-spins. Least-square fits to the data indeed indicate a drop in the B(E2) values 
towards the high-spin end. However, the error bars are too large to reach an unambiguous 
conclusion whether this effect is genuine or not. 

V. CONCLUSIONS 

In this paper, we have reviewed a microscopic basis and a practical formulation of the 
IBM for application to superdeformed nuclei. The availability of analytical formulas owing 
to the 1/N expansion technique means fast and efficient analysis of data. As first examples, 
we have considered the superdeformed bands in the Hg-Pb and Gd-Dy regions. A good 
description of data is obtained in the Hg-Pb region confirming the simple quadrupole nature 
of these superdeformed bands. In the Gd-Dy region, the dynamic moment of inertia exhibits 
large variations which can not be accommodated in a simple collective model. Such varia- 
tions are due to the single particle degree of freedom and require extension of the model for 
a better description of the data. Finally, the formalism can be used in investigating some 
other interesting features of superdeformed nuclei such as identical bands and C4 symmetry 
which will be pursued in future work. 
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FIGURES 

FIG. 1. Occupation probability of each spherical basis in the Nilsson potential: (a) protons and 
(b) neutrons of 194 Hg and (c) protons and (d) neutrons of 152 Dy. Two values of the deformation 
parameter are considered, S = —0.13 (dashed line) and 0.40 (solid line) for 194 Hg, and 5 = 0.25 
(dashed line) and 0.50 (solid line) for 152 Dy. The spherical magic numbers are indicated. 

FIG. 2. Comparison of the 1/N expansion results (lines) with the exact numerical ones (circles) 
for N = 10 bosons. Fig. 2-a shows the dynamic moments of inertia obtained from the third order 
calculation with VAP (solid line), third order with VBP (dashed line), and second order with VAP 
(dotted line). The parameters of H are k = 20 keV and q = 0.7. Fig. 2-b compares the E2 matrix 
elements. 

FIG. 3. Systematic behaviour of dynamic moment of inertia for various values of the 
quadrupole parameters c/22> 924 and §44. In each figure, two of these three parameters are fixed at 
the SU(3) value while the other takes 0.6, 0.8, 1.0, 1.2, 1.4 times the SU(3) value. N = 30 and n 
is chosen such that the moment of inertia are all normalized to 100/i 2 MeV~ 1 at u = 0. 

FIG. 4. Effect of the simultaneous scaling of the quadrupole parameters (a) on dynamic moment 
of inertia, and (b) on B(E2) values. The dynamic moment of inertia curves are normalized to 100 
ft 2 MeV _1 at oj = and the B(E2) values are normalized to B(E2; 2 -> 0). 

FIG. 5. Comparison of the experimental dynamic moment of inertia in 190 ~ 194 Hg (circles) with 
the super IBM calculations (solid lines). N = 29, 30, 31, « = 35, 34, 33 keV, q = 0.68, 0.72, 0.72 are 
taken for 190-194 Hg, respectively. The data are from [pTof] . 

FIG. 6. Same as Fig. 5 but for 192 ~ 19 6pb. N = 30, 31, 32, k = 33, 33, 34 keV, q = 0.66, 0.65, 0.67 
are taken for 192 ~ 196 p D) respectively. 

FIG. 7. Same as Fig. 5 but for 148 - 150 Gd. N = 29, 30, k = 41, 27 keV, q = 1.35, 1.73 are taken 
for 148 ~ 150 Gd, respectively. 
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FIG. 8. Same as Fig. 5 but for 152 - 154 Dy. N = 31,32, « = 42,43 keV, q = 1.07, 1.02 are taken 
for 152_154 Dy, respectively. 

FIG. 9. Comparison of the experimental B(E2) values in 192 - 194 Hg (circles) with the super IBM 
calculations (solid lines). The boson effective charges are e = 0.140 and 0.124 eb for 192 ^ 194 Hg, 



respectively. The data are from [15,22]. 
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TABLES 

TABLE I. Probability (%) of each angular momentum component in the yt-pair. The /1-pair 
is obtained from the Nilsson + particle number conserving BCS wave function. Two cases of 
5 = —0.13 and 0.40 are shown for 194 Hg. The probability of each boson in the intrinsic boson of 
the IBM in the SU(3) limit is also listed for comparison. 
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